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Thermal (or pseudo-thermal) radiation has been recently 
proposed for imaging and interference types of experiments to 
simulate entangled states. We report an experimental study 
on the momentum correlation properties of a pseudo-thermal 
field in the photon counting regime. The characterization 
and understanding of such a light source in the context of 
two photon physics, especially its similarities and differences 
compared to entangled two-photon states, is useful in gaining 
knowledge of entanglement and may help us to assess the real 
potential of applications of chaotic light in this context. 



The study of thermal radiation has been historically 
important in optics [1]. Before the invention of the laser, 
thermal sources were probably the only light sources 
available in the laboratory for optical experiments and 
observations. Although coherent light has taken the ma- 
jor role in modern optics laboratory, thermal, or pseudo- 
thermal, sources are still important for the understanding 
of fundamental concepts, such as optical coherence and 
photon statistics [1-4]. 

Recently it has been theoretically proposed to use ther- 
mal, or quasi-thermal, radiation [5] in the context of two- 
photon imaging and interference [6]. The proposal fits 
in the lively ongoing discussion about the possibility of 
using classically correlated light sources to simulate the 
behavior of entangled photon states. In the past two 
years, in fact, successful attempts to simulate some of 
the surprising features of entangled states with classi- 
cally correlated light have been made [7], leading to an 
interesting discussion about which features are peculiar 
to entanglement, and to what extent a classical corre- 
lation can simulate them [8]. Furthermore it has been 
argued that thermal light exhibits exactly the kind of 
correlations necessary to simulate, at least qualitatively, 
most of the results obtained in the past through entan- 
gled photon states [5]. 

In this paper we report an experimental study on the 
momentum correlation properties of thermal radiation 
in the photon counting regime. A theoretical analysis is 
also proposed to clarify the physics behind the observed 
phenomenon, especially the similarities and differences 
with entangled two-photon states. It may not be quite 
correct to use the terminology "correlation" to describe 
the physics behind the experimental observation. In fact, 
the observed result comes naturally from the complete 
absence of correlation, as is hinted by the same name of 
chaotic light. 

In principle the term "thermal radiation" should refer 



only to radiation coming from a blackbody in thermal 
equilibrium at some temperature T. But in reality, some 
characteristics of true thermal fields, like the extreme 
shortness of their coherence time, have imposed serious 
obstacles to their use in actual experiments, and there- 
fore since the early days of quantum optics there has 
been a great interest in the realization of more utilizable 
sources that could simulate the behavior of true thermal 
fields (gas discharge lamps, randomized lasers etc.). We 
usually describe these kind of sources as pseudo-thermal; 
they are actually all chaotic light sources that can be 
modelled as a collection of independent atoms emitting 
radiation randomly [9] . The most commonly used among 
the pseudo-thermal fields is the one developed by Mar- 
tienssen and Spillcr [10] in the mid of 1960 's. The prin- 
ciple of the generation is very simple: coherent laser ra- 
diation is focused on a rotating ground glass disk so the 
scattered radiation is chaotic with a Gaussian spectrum. 

We start by describing the light source used in this ex- 
periment. We experimentally studied the characteristics 
of the source to be sure that it indeed behaved like a 
thermal field. To this purpose we used an experimental 
setup similar to the one that was later used to perform 
the momentum correlation measurement (Fig. 1). First 
we focused on the temporal behavior and repeated a his- 
torical experiment of Arecchi et al [11]. The point of 
the experiment was to study the time distribution of the 
coincidence counts coming from the source under study. 
Fig. 2 shows the results of the measurements. Basically 
we plot a histogram of the number of coincidence counts 
versus the time difference of the photon registration clicks 
from the two detectors. The temporal behavior of our 
source is similar to the one obtained with light coming 
from stars and gas discharge lamps [12-14], but with a 
higher coherence time (around 1/is, mainly determined 
by the rotation speed of the disk and the properties of 
the grit of the ground glass). In practice this graph is also 
useful for "calibrating" the coincidence time window; in 
fact in the rest of the experiment we set a window of 
about 600 ns around the peak of Fig. 2 and measured 
the number of coincidence counts in that window. Then 
we studied the spatial behavior of the coincidence counts 
using the setup in Fig. 1. The point of this measurement 
is to study the behavior of the spatial correlations when 
the size of the source is varied. This kind of measure- 
ment is basically the one used to determine the angular 
size of the stars in Hanbury-Brown and Twiss stellar in- 
tcrferomctry. Experimentally we fixed detector D\ in one 
position and scanned horizontally the position of detector 
Di for different sizes of the source obtained by focusing 
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or defocusing the Hc-Ne laser onto the rotating disk. The 
results are shown in Fig. 3. As expected, the bigger is 
the angular size of the source, the smaller is the width of 
the correlation function, in agreement with stellar inter- 
ferometric predictions. 

These two measurements confirmed that our source 
simulated both spatially and temporally the behavior 
of thermal radiation. We then proceeded on measuring 
the momentum-momentum correlation of the field. The 
setup for the measurement is shown in Fig. 1. After the 
source, a lens (L) of focal length / = 75cm was inserted, 
and then the radiation was split by a non-polarizing cube 
beam splitter {I = 25mm). In the focal plane of the lens 
L we put the detection system composed of two horizon- 
tally displaceable fibers (60^m diameter) connected to 
two single-photon counting modules. The output pulses 
from the two detectors were then sent to an electronic 
coincidence circuit. The setup is designed to study the 
momentum properties of the fields, in fact the lens L 
was close to the source and the detectors were placed 
in its focal plane to ensure the far-field condition. In 
this situation, in principle, the fiber tip (detector) scans 
the transverse wavevector plane, or transverse momen- 
tum plane, of the field. Therefore by scanning the fibers 
in the horizontal direction continuous information about 
the behavior of each k-vector (momentum) of the field 
will be obtained. 

We performed three sets of measurements for three dif- 
ferent position of detector D\ ; in every measurement we 
kept the position of D\ fixed and scanned horizontally 
the position of detector D 2 . The results are shown in 
Fig. 4. It is clear that the correlation function "followed" 
the movement of detector D\\ a shift in the position of 
D\, causes a shift in the same direction of the correlation 
function. The natural deduction from this measurement 
is that the field is indeed correlated in momentum: in 
fact every time a position of Z?i is fixed, one k-vector 
from the source is selected, therefore the behavior of the 
coincidence counts show that only one k-vector arriving 
at D 2 is correlated to the selected k-vector arriving at 
D\. In the central part of Fig. 4 the single counts are 
reported, notice that (1) the single counts are flat over 
the entire range of detection suggesting the completely 
random transverse distribution of the field; (2) the low 
level of photon counts (60000 per second) suggests that 
we were effectively operating in the condition of only two 
photons present in the time duration of interest. 

The phenomenological conclusion of these measure- 
ments would then be that the pseudo-thermal field we 
have measured is indeed correlated in momentum. This 
is a property that makes this source particularly interest- 
ing to study in comparison with entangled photon pairs 
produced in Spontaneous Parametric Down Conversion 
(SPDC) in the context of two-photon imaging and inter- 
ference features [6,5]. However, the predicted, and now 
experimentally measured, analogy between the effects ob- 
tained with entangled photons and thermal radiation has 
to be analyzed very carefully. In order to clarify the point 



let us justify theoretically the measured experimental ef- 
fect. The quantity that rules the experiments involving 
coincidence measurements is the second-order Glauber 
correlation function [15]: 

G< 2 >(ti,x i; t 2 ,x 2 ) = (1) 
= Tr[pE { -\t 1}Xl )E { 2 -\t 2 , X2 ) 

Where is a vector in the plane of the i th detector 
perpendicular to the direction of propagation of the ra- 
diation, p represents the density matrix of the quantum 
state under consideration, and E^iti,^), i = 1,2, are 
positive-frequency and negative- frequency components of 
the field at detectors D\ and D 2 that can be written as: 

E\ + \U^) = J2 y^ 9i (x,,k)a k e-^. (2) 
k 

Here ^(x^k) stands for the Green's function associated 
with the propagation of the field towards the i th detec- 
tor [16]. Plugging the expression of the electric fields in 
Eq. 1, we obtain the equal time second-order Glauber 
correlation function in the following form: 

G (2) (x i; x 2 )oc 

EEEE^( x i' k )ff2(x 2 ,k') 

k k' k" k " 

x g 2 (x 2 , k")ffi (xi , k"')Tr [pa^a^ a k " a k >» ] 

Let us examine the last term of Eq. 3, in particular the 
density matrix of the state of the radiation we are study- 
ing. Usually the density matrix that is used for single 
mode thermal light is expressed in the number state ba- 
sis as follows [9] [17]: 

p = V - \n)(n\ (4) 

' ^ (1 + n ) 1+ ™ n 1 w 

In the multimode case we can take {n k }} as basis vec- 
tors, where the symbol {n- k } indicates a set of number 
of photons n ki' n k 2 ' n k 3 e * c - m ^ e m °dcs ki,k 2 ,k 3 re- 
spectively, and {n k }} are product states of the number 
states n^} of each mode k. In thermal and chaotic light 
the different modes of the field are independent [9] , thus 
the total density matrix will be the product of the con- 
tributions from the different modes as follows: 

P = E l{»k»({»k}l II 77T7 k a + » k (5) 

where (n k ) represent the average number of photons in 
mode k. 

A comparison between the density matrix that de- 
scribes the thermal radiation and the state of the entan- 
gled photons from SPDC can readily clarify the intrinsic 
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difference between the two sources. While in the ther- 
mal case it is intuitive and apparent from Eq. 5 that the 
different modes of the radiation are completely indepen- 
dent, it is immediately clear from the state of the SPDC 
radiation that there is perfect correlation in momentum 
exhibited by the source. 

The trace in Eq. 3 turns out to be: 

T^a^a^a^] = <n) 2 (5 kk » V,k'" + 5 k,k'" 5 k',k") 

(6) 

where we have assumed the average number of photons 
in each mode to be constant (as is confirmed by the single 
counts behavior in Fig. 4). Plugging the previous result 
in Eq. 3 we obtain the following expression for the second- 
order correlation function: 

G^(xi;x 2 )oc (7) 
J AlSifxj.q) | 2 |rf 2 q|.9 2 (x 2 ,q) | 2 

+ I J d 2 qffi(xi,q).92(x2,q) | 2 • 

Where q is the transverse component of the k-vector. 
Notice that we have converted the summation of Eq. 2 
into an integral . It is easy to see that J ePq | <7i(xj, q) | 2 
will actually lead to the number of single counts n, of the 
i th detector. The expression can be further simplified: 



<a k °k'> = (n)5 kk , 



(11) 



G(2)(xi;X2) K ! + iKqrifo.qtofo.q) 



nin 2 



(8) 



Inserting the Green's functions associated to the optical 
setup into Eq. 8 we obtain: 



G< 2 >(xi;x 2 ) cx 1 + 



■ Ji(7ra|xi - Xgj/A/) 
- 7 ra|x 1 -x 2 |/2A/ 1 



(9) 



where J\ is the Bessel function of the first order and a 
is the size of the source. This result explains clearly the 
experimental observations of both Fig. 3 and Fig. 4. In 
fact, it predicts that the bigger the size of the source, 
the smaller is the width of the correlation, and above 
all it predicts that shifting the position of detector D\ 
will cause the observed shift of the correlation function. 
Notice that, due to the first term in Eq. 9, the maximum 
achievable visibility in these kind of experiments is 33%. 

Although the theoretical calculations are able to pre- 
dict the experimental observations, the physical inter- 
pretation of the experimental results and mathematical 
formalism remains yet to be fully explained and under- 
stood. To clarify this point, let us focus for a moment on 
the following quantities: 

(a\a k ,)=Tr[pa\a k ,\ (10) 

( a k a k' a k" a k"') = T K/^a k /a k »a k <"] 

Substituting into this equation the expression for the den- 
sity matrix that describes our source we obtain: 



H a k' a k" a k"') = ( n ) ( 5 k,k" 5 k',k"' + (5 k,k"" 5 k\k") ( 12 ) 

Notice that relations of the kind in Eq. 11 are gener- 
ally used in statistics to express the absolute statistical 
independence of the variates under study. In this case 
therefore Eq. 11 can be interpreted as the mathemati- 
cal expression of the lack of correlation between different 
modes of the field [1] and the experimental observations 
would be a natural manifestation of this property. 

As far as Eq. 12 is concerned, in the context of entan- 
gled photon states, a formally similar equation is usually 
considered explicative of the perfect momentum correla- 
tion between the two entangled photons due to the phase 
matching conditions of the nonlinear process inside the 
crystal. The interpretation of this mathematical relation 
in terms of correlation between k-vectors could be mis- 
leading in the case of chaotic light. In fact, it would 
employ a terminology peculiar to entangled photons for 
a process that is intrinsically different. 

Besides the physical interpretation, however, the the- 
oretical and experimental results reveal that the source 
under study, and thus thermal radiation in general, is 
promising for two-photon interference type experiments, 
except for the lower visibility (33% maximum) [18]. From 
a fundamental point of view it can represent a perfect tool 
for studying the differences and/or the analogies between 
quantum entanglement and its classical simulation. 

In conclusion, we have experimentally studied the far- 
field momentum "correlation" of chaotic light fields and 
we have provided a theoretical explanation of the results 
with a physical interpretation of the phenomenon. The 
observed correlation is in fact the experimental evidence 
of the complete absence of correlation between different 
modes of the radiation. The experiment can help us to 
assess to what extent thermal or pseudo-thermal fields 
might be employed to simulate certain features of entan- 
gled photons in the context of two-photon imaging and 
interference. 
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FIG. 1. Sketch of the experimental setup 



FIG. 2. Histogram of number of joint detection counts vs 
time difference of the two photo-electron events. The size of 
each channel is 0.3 ns. 



FIG. 3. Normalized second order correlation function vs 
position of detector D2. The various sets of data (points) 
are for different sizes of the pseudo-thermal source (circles 
~ 150^im; squares ~ 350/im; triangles ~ 550pm; diamonds 
~ 1100/im) and agree with the theoretical expectation (solid 
lines). 



FIG. 4. Normalized second order correlation function vs 
position of detector D2 for three fixed positions of detector 
Di: (a) xi ~ —1.75mm; (b) x\ ~ 0mm; (c) x\ ~ 1.75mm 
(source size ~ 550/im). In part (b) is shown also a plot of 
single counts vs position of the detector for both D\ (filled 
squares) and D2 (hollow squares). The level of the counts is 
around 60000 per second. 



Figure 1. Giuliano Scarcelli, Alejandra Valencia, and 
Yanhua Shih. 
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Figure 2. Giuliano Scarcelli, Alejandra Valencia, and 
Yanhua Shih. 
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Figure 3. Giuliano Scarcelli, Alejandra Valencia, and 
Yanhua Shih. 
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Figure 4. Giuliano Scarcelli, Alejandra Valencia, and 
Yanhua Shih. 
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